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Abstract 

The duality transformation is applied to the Fisher zeroes near the ferromagnetic critical 
point in the q > 4 state two dimensional Potts model. A requirement that the locus of the 
duals of the zeroes be identical to the dual of the locus of zeroes in the thermodynamic limit 
(i) recovers the ratio of specific heat to internal energy discontinuity at criticality and the 
relationships between the discontinuities of higher cumulants and (ii) identifies duality with 
complex conjugation. Conjecturing that all zeroes governing ferromagnetic singular behaviour 
satisfy the latter requirement gives the full locus of such Fisher zeroes to be a circle. This locus, 
together with the density of zeroes is then shown to be sufficient to recover the singular part of 
the thermodynamic functions in the thermodynamic limit. 



^Supported by EU TMR Project No. ERBFMBI-CT96-1757 
^ Current address 



Introduction The q-state Potts model ]^], introduced in 1952 as a generalization of the Ising 
model 1^] , has become the generic model for the analytical and numerical study of first and second 
order phase transitions . Apart from the one dimensional case [j^] , the only solution which exists 
to date is for the q = 2 (Ising) model in two dimensions and in the absence of an external magnetic 
field IQ]. The partition function for the standard Potts model is Zl{(3) = J2{ai} ^^P {P'l2(ij) ^o-io-j) 
where P = l/lksT), and T is the temperature. The spin cjj at site i on a d-dimensional lattice takes 
values 1,2, ... ,q and the total number of sites is V = L'^. Despite the absence of a full solution for 
general q, some exact and, indeed, rigorous results are obtainable in d = 2. The first of these is 
that, up to an irrelevant multiplicative constant, the form of the partition function is unchanged 
under a duality transformation in the thermodynamic limit In terms of the (low temperature 
expansion) variable u = e~l^ , this duality transformation is u ^ T^{u), where 

The critical temperature at which the phase transition occurs is invariant under (0) and is given 

by i, 

1 

V(uc) = Uc or Uc = — ■ (2) 

\ cj c ^ W 

Baxter has shown that at the critical point the model is equivalent to a solvable homogeneous ice- 
type model js], ^ . By deriving the latent heat at criticality it was shown that the phase transition in 
the two-dimensional model is first (second) order for q > 4: (q < 4). In fact, for the q > 4 case, the 
exact values of the latent heat, the mean internal energy and the specific heat discontinuity (but 
not, for example, the mean specific heat) are known [|^, ||, ^, ^. The full form of the free energy 
(and derivable thermodynamic functions) of the Potts model has, however, never been calculated for 
general q and general T. In the words of Baxter, solving the Potts model for general temperatures 
is, therefore, "a very tantalising problem" |Q. The Potts model is reviewed in 

In this paper the problem is approached using a remarkably general and recently derived result 
concerning the partition function zeroes of models with a first order phase transition ||8| . 

For finite systems the zeroes of the partition function Q are strictly complex (non-real). As L ^ oo 
one expects these zeroes to condense onto a smooth curve whose impact on to the real parameter 



axis precipitates the phase transition. Fisher [10| emphasized the application of zeroes in the 
complex temperature plane to the study of temperature driven phase transitions. In particular, 
in Ijl^, the Kaufman solution ||ll| of the two dimensional Ising model was used to show that the 
Fisher zeroes (also called complex temperature zeroes [jl2|) are dense on two circles in the complex 
u-plane in the thermodynamic limit. 

Based on similarities with the Ising case, Martin and Maillard and Rammal |13, 14] conjectured 
that the locus of Fisher zeroes in the d = 2, g-state Potts model be given by an extension of (^) 
to the complex plane, namely I'(n) = u* where u* is the complex conjugate of u (although, on 
this basis alone 'it is not clear where this requirement comes from' |]13[). This identification yields 
a circle with centre —l/{q — 1) and radius y/q/{q — 1). When q = 2 this recovers the so-called 
ferromagnetic Fisher circle of the Ising model [llO||. In the Ising case, the partition function is 



actually a function of n^. There, the second (so-called antiferromagnetic) Fisher circle comes from 
the map u — > (/? — > — /?). Numerical investigations for small lattices at g = 3 and 4 [l^, 0] 
provided evidence that the Fisher zeroes indeed lie on the circle given by the identification of duality 
with complex conjugation. However the numerics are highly sensitive to the boundary conditions 
used and the situation far from criticality remained unclear. Some progress was recently made in 
the non-critical region using low temperature expansions for 3 < g < 8 ||12|. 



Recently, and on the basis of numerical results on small lattices (up to L = 7) with q < 10, it has 
again been conjectured that for finite lattices with self-dual boundary conditions, and for other 
boundary conditions in the thermodynamic limit, the zeroes in the ferromagnetic regime are on 
the above circle The conjecture of |15| was, in fact, proven for infinite q in [^]. This circle- 
conjecture is similar to another recent conjecture |17], namely that the Fisher zeroes for the (/-state 
Potts model on a triangular lattice with pure three-site interaction in the thermodynamic limit 
(which is also self-dual ||l^) lie on a circle and a segment of the negative real axis. 

All of the above conjectures regarding the locus of Fisher zeroes rely, at least in part, on numerical 
approaches. In this paper, the problem is addressed analytically. A requirement that taking 
the thermodynamic limit and application of the duality transformation to the Fisher zeroes be 
commutative in the q > 4 case (i) recovers the ratio of specific heat discontinuity to latent heat and 
corresponding relationships between the discontinuities of higher cumulants and (ii) analytically 
identifies duality with complex conjugation. Conjecturing that all zeroes governing ferromagnetic 
singular behaviour satisfy the latter requirement, the locus of such Fisher zeroes is shown to indeed 
be a circle. This locus, together with the density of zeroes is then shown to be sufficient to recover 
the singular form of all thermodynamic functions in the thermodynamic limit. 



Thermodynamic Functions Consider, firstly, the finite-size system. For finite L, the partition 
function can be written as a polynomial of finite degree in u, and as such, can be expressed in terms 
of its complex Fisher zeroes Uj{L) as Zl{P) oc 0^=1 (u ~ Uj{L)). The free energy is defined by 
PfiP) = — l'nZ{P)/V. The internal energy is therefore 

The specific heat and the general n}^ cumulant are respectively defined as 
Using the notation 

P/ Pc P \Pc 

for the discontinuity in the n^^ cumulant at the critical temperature, the exact results |5|, ^, |^] 
(in the thermodynamic limit) are 



Ae 
Ac 



2(1 + -^"j tanh Y[ tanh^ (nO) , 



n=l 



(7) 
(8) 



where = In y^Jq/A + ^q/^ — Ij . Further results include the general higher cumulant combination 
ln{Pc) - {-Tln{f3t) determinable from duality U, M. 



Partition Function Zeroes Recently, Lee |^] has derived a general theorem for first order phase 
transitions in which the partition function zeroes can be expressed in terms of the discontinuities 
in the thermodynamic functions (for finite size as well as in the infinite volume limit). For a system 
with a temperature-driven phase transition, Lee's result for the Fisher zeroes is 

Ing , , (2j-l)7r _ , ^ (A^ A7n 

VAe FAe " '^^^ + 2! Ae + n\ Ae 

where the reduced temperature is t = 1 — /3//3c- Inverting, we find ||8| 

/3cRetj(L) = + ^3/4 + + . . . + , 

±/3cImtj(L) = + ^2/| + ^4^1 + . • • + , (10) 

where Ij = (2j — l)7r/(yAe) and 0{\/V) represents terms which vanish in the infinite volume limit 
and where the coefficients An are easily calculable, the first few being [H] 



A, = , (11) 



A, . _5.3.5..^-^ . ,13) 
A, = 144-21^;|?i+3f^)%6A.^-^ . (14) 



3!Ae V3!Ae/ 4!Ae 5!Ae 

A = 424-84^?4?^ + 28Aif4^V + 28A?4ji-74^4ji-7yli4^ + 4^(15) 
^ ^3!Ae V3!Ae/ M!Ae 3!Ae 4!Ae 5!Ae 6!Ae ^ ^ 



The Locus of Zeroes From (10) the real part of the zeroes (in the thermodynamic limit) can 



be expressed in terms of their imaginary parts as (3cRet = £(/3cImt) where 

C{e) = Aie^ + {-2AiA2 + A3)0^ + i7AiAl-2AiAi-4A2A3 + A5)e^ + ... . (16) 

The zeroes are thus seen to lie on a curve. In the complex u upper half-plane the equation of this 
curve is 

^(+)(e)=^,e^W+^^ . (17) 
This defines the locus of zeroes in the infinite volume limit. 



The Dual of the Locus of Zeroes Applying the duality transformation (Q) to 7^'^^(^) and 
expanding in powers of 9 gives 



ReV (-f'^+\e) 



-72Aiq + 72Aiqi - l2Aiq^ + 24^?gi - l2Alq^ - 48^1^2?^ + 2AA-iq^) + 



(18) 



Imp (7^+^(6') 



6g 



6 - 6g2 + g - 12^1^2 + QAiq 



120g2 



-120 + 240g^ 



-150g + 30^2 -q^ + 480ylig2 - 720^ig + 280^1^2 - 20Aig^ - 360^?g + 360A?g2 



-mA\q^ + 480^1 Asgt - 240^1^2^2 - 240^3gt + 120^3^^) 



+ ... 



(19) 



The Locus of the Duals of the Zeroes Alternatively, applying the duality transformation (|i|) 
directly to the j*^ zero in the finite-size system and expanding again, gives 



/3cRetf(L) 



Afl] + Af^lf + Afl^ + ... + 0{l/V) , 



TPcluityiL) = Ij + A^lj + A^I^ + ... + 0{l/V) 
where the first few coefficients A^ are 



(20) 



A? 
A^ 



^4 



^5 



-1 



+ 2q 2A1 + A2 



Ai , (21) 

(22) 

+ ^~^Ai-q-^A\ + 2q-^2A2- A^ (23) 
-q-^ + - \q'^Ai - Aq-^Ai + ^"^Al - 3g~^yl2 + 2g-^^i^2 + 2g"^^3 + ^4 (24) 



q 2 

12^ 



_ 3 
q 2 



3 

1 _1 1 _3 

q 2 + -q 2 + 

360 2^ 



t - ^q'^Ai - bq-'^Ai + Qq'l Al + Al - q-^A{ - ^9^^^ 



-4g 5^2 + 6g ^^1^2 + 9 ^^l + Sg ^^3 - 2g 5^1^3 + 2^ -A5 



2 

(25) 



Prom ( pO| ) the real part of the dual zeroes can be expressed in terms of their imaginary parts in 
the thermodynamic limit as /JcRetj^ = C^ilmt^) where 

C^{e)=A^e'^ + {-2A^A^ + A^)e^ + {lA^Af -2AfA'i -AA^A^ -A^)e^ + ... . (26) 

Therefore, the locus of the dual of the upper half-plane zeroes in the thermodynamic limit is given 

by 



^(+)^(0)=^,e^"W-^ 
The expansion of this locus of duals is 

Re7( 



(27) 



v(+)''(0) 



1 + 



^ ' 1 + 2^f ) + ^ (1 - 12^f - A8AfA^ + 2AAf + 12^f ^) + 



-( 

2! V 



Im7(+)^(6') 



D , 



-Ur 



II (-1 + 6^f ) + ^ (1 - 20^f + 60Af ^ - 240^f + 120^^ ) + 



\D aD 



(28) 
.(29) 



Identification of the Dual of the Locus with the Locus of the duals In deriving the dual of 
the locus of zeroes (^) and (^), the duality transformation was applied after the thermodynamic 
limit of the positions of the zeroes (i.e., their locus) was taken. In ( pSD and ( p9| ) the duality 
transformation was applied to the zeroes before taking the thermodynamic limit. Even in the 
case where the finite-L system does not have duality-preserving boundary conditions, taking the 
thermodynamic limit restores self-duality. The dual of the (thermodynamic limit) locus of zeroes 
and the (thermodynamic limit) locus of the duals of the zeroes must be identical. We demand, 
therefore, that 

V = (30) 

order by order in the expansion in 9. Up to 0{9^) this is trivial. To O{0^) and (separately at) 
0{d'^) they are identical if Ai = 1/(2^). From (0), this is the result (||). The identity (||) at 
0{9^) and (separately at) 0{e^) gives A3 = A^l ^ - q-'^/'^{q - 3)/24, which from ^ and (|l|) 
means that 

A74 = Aa73 + i-^Ae . (31) 

Higher order results are obtainable using a computer algebra system such as Maple. To orders 
0{e'^) and 0{e^) and (separately) to orders 0(0^) and 0{e^^) one finds 

A76 5 A75 g-20 A73 1 1 ^ 1 



6!Ae 2gi/2 5!Ae 8g3/2 3!Ae 6!gi/2 8g3/2 2g5/2 



and 



A78 ^ 7 A77 / 5 35 \ A75 / 3 15 21 \ A73 

8!Ae 2gV2 71 Ae \2Aq^ /2 4g3/2 ) 51 Ae Ulg^/^ 16g3/2 + 2g5/2 j 3!Ae 

1 23 5 17 

^8V72 ~ 960^ ^ 8^ ~ 8^ ' ^ ' 

respectively. These results and further results for the higher cumulants at criticality are also 
obtainable directly from the duality transformation (|3l|) (see 19, [20[]). 



The Full Ferromagnetic Locus of Zeroes Putting the above equations into (pi[) - (p5[) (and 
their higher order equivalents) yields A^ = Aj (this has been verified up to j = 8). Therefore (at 
least up to 0^^) the dual of the locus of zeroes is the complex conjugate of the original locus of 
zeroes. We now assume that this is the case for all 9. Then, the full ferromagnetic locus of zeroes 
(that part of the full locus which intersects the real temperature axis at the physical ferromagnetic 



critical point) is found by identifying |13, 14] 

p(7(+)(0) 



7 



i+r 



(0) 



(34) 



where 7^+)* represents the complex conjugate of 7^+). The full ferromagnetic locus is then |T^, 0] 



7(0) 



-A- 



■1 + Vq^ 



id 



(35) 



This circular locus is analogous to the circle theorem of Lee and Yang |^. In the field driven case, 
where one is interested in Lee- Yang zeroes in the complex z = exp h plane {h is an external magnetic 



field), formulae analogous to (^) - ( [l5| ) apply where e, c, etc. are replaced by the corresponding 
derivative of the free energy with respect to h (the magnetization m, the susceptibility X: etc.). 
There, the partition function is unchanged under h — > —h and consequently A^i = for even 
I. Therfore C{9) = and the locus of zeroes is z = expiO. This is Lee's proof of the Lee- Yang 
theorem One observes that considering h ^ —h as a self-duality map and identifying it with 
complex conjugation yields this locus. 



The Singular Parts of the Thermodynamic Functions in the Thermodynamic Limit 

Prom (^), the density of zeroes in the temperature driven case is [^] 

(In {{Vq + 1)7(^)))' + (1^ ((v^ + 1)7W))' + • • •} • (36) 

The internal energy is (from (^) or |2^, |22[| ) 

e = cnst. + u ^^^\^, de . (37) 

Therefore, from (|35|), ( p^ ) and (^), the internal energy is 

e(/3</3c) = eo (38) 

e(/3 > /?.) = eo - Ae |l + ^^(/3. - /?) + i^(/3. - /5)^ + i^(/3. - /?)3 + . . .} (3<)) 

with Co a constant (one expects that when separate Fisher loci which don't cross the positive real 
temperature axis are accounted for, eo becomes temperature dependent). At fic the internal energy 
discontinuity e(/3 = /3~) — e{(3 = /?+) = Ae is recovered. Appropriate differentiation recovers the 
discontinuities in specific heat and higher cumulants. Thus the full locus (^) and the density ( |3^ ) 
are sufficient to give the singular parts of the thermodynamic functions in the infinite volume limit. 



Conclusions In summary, we have applied the duality transformation (|T|) , under which the d = 2 
g-state Potts model is invariant, to the Fisher zeroes recently found by Lee for systems with a 
first order phase transition. The requirement that the dual of the locus of zeroes be identical to the 
locus of the duals of the zeroes in the thermodynamic limit (i) recovers the ratio of specific heat to 
internal energy discontinuity at criticality and the relations between the discontinuities of higher 
cumulants and (ii) identifies duality with complex conjugation. 

Conjecturing that all zeroes governing ferromagnetic critical behaviour satisfy (ii) gives that this 
full locus is the circle (35) in the complex u-plane. The equation ( |35| ) was first conjectured by 
Martin and Maillard and Rammal |]l^, |lj] on the basis of analogous Ising results [10|. The same 
conjecture, based on numerical results for small lattices was made in |15] and proven for infinite q 
in [^] (see also ^ 



The locus (|35| ) , together with the density of zeroes is sufficient to recover the singular parts of all 
thermodynamic functions in the thermodynamic limit. It is to be expected that the regular parts 
come from separate loci of zeroes which don't cross the positive real temperature axis. 
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